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Some Remarks on Lambert's Theorem . By H. C. Plummer, M.A. 

1. It is a curious fact that the proofs of a theorem so celebrated 
as Lambert’s, as they appear in at least the majority of text-books, 
are very unsatisfactory in the treatment given to the ambiguities of 
sign which must occur in the statement of the theorem. In most 
of the cases which arise in practice the arc of the orbit is com¬ 
paratively short, and there is really no question about the signs to 
be chosen. This, coupled with the fact that it was commonly 
supposed that the form of the theorem was in general unsuitable 
for use in the determination of orbits, has apparently caused a 
certain carelessness in the demonstration. Recent work has, how¬ 
ever, tended to place the value of the theorem more clearly in 
evidence.* For the case of elliptic motion the essence of the matter 
in regard to the distinction in the ambiguous signs has been' given 
by Callandreau. f For hyperbolic motion under attraction to a 
centre a simple proof was given by Adams J in terms of hyperbolic 
functions, but the question of signs was ignored by him altogether. 
The case of motion under a repulsive force along the convex branch 
of an hyperbola does not seem to have been discussed at all, though 
the motion of matter in a comet’s tail may give this case a real 
practical importance. It may therefore be useful to give here some 
examination to the methods of proof in the several cases, with such 
remarks as may be suggested incidentally. 

2. The case of elliptic motion will be discussed first. Let 
E 1? E 2 be the eccentric anomalies of the starting and end points 
respectively, and let E 2 >E X , both angles being less than 360°. In 
the first instance we suppose that the motion is in the same sense 
as the anomaly increases. It is easy to prove that, yu, being the 
force at unit distance, 

^-a~U = e - 8 - (sin € - sin 8 ) 

where 

sin 2 = (r x + r 2 + c)j^a = m x 2 , sin 2 = (r a 4- r 2 - c)!^a =m 2 2 

The relations of c and 8 to the eccentric anomalies are given by 
e - 8 = E 2 - E x , cos -1- 8 ) = e cos J(E 2 + E x ) 

The angle |(c + 8 ) is merely defined by the latter relation, and 
may be taken between the limits o° and 180°; and £(* - 8 ) must 
lie between the same limits. Hence e must lie between o° and 
360°, while 8 must lie between ± 180°. It follows that sin Je and 
cos JS are both positive. If, then, e x and are the principal values 
resulting from 

sin Je = 4- m 1 , sin = + m 2 

* Cf. M.N., lxiii. p. 147. 

t Ammles de I’Obs. de Paris, Memoires, vol. xxiii. 

X Collected Scientific Papers , vol. i. p. 411. 
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we have to choose between e l and 2ir - e x and between ± 8 r Now 
the intercept x on the major axis between the chord joining the 
end points and the centre is given by 

x = a sec |(E 2 + E x ) cos ^(E 2 - E x ) 

or 

(x + ae)/(x - a) = {cos £(E 2 - E x ) + e cos J(E 2 + E^} 

-=■ {cos KE 2 - Ej) - cos |(E 2 + Ej)} 

= cos |e cos £ 8 / sin JE 1 sin ^E 2 

Under the conditions imposed, the sign of cos Jc is the same as 
that of the fraction on the left: this shows that e = e 1 or 27r - c x 
according as the chord intersects the axis outside or inside the 
range between the second focus and “ perihelion.” Similarly 

(x - ae)/(x - a) = sin Je sin £ 8 / sin £E 2 sin JE X 

which shows that 8 = S 1 or - 8 X according as the chord does not, or 
does, intersect the axis between the attracting focus and 4 ‘perihelion.” 

3. Let now 

e x - 8 X - (sin e L - sin 8 : ) = T x = nT 1 
€1 + Si - (sin € 1 + sin 8jj = /Jet ~f T 2 = nT 2 

Then, under the conditions assumed, we have found that 

t = T 2 , 27 rjn - T 2 or 2 tt In - T 2 

according as the segment of which the curved boundary is described 
in the motion contains (1) no focus, (2) the empty but not the 
attracting focus, (3) both foci. But the conditions have precluded 
motion past “perihelion,” and this condition must be removed. It 
is merely necessary to consider the motion over the complementary 
parts of the ellipse corresponding to the cases already distinguished. 
We thus find that 

t= 27 r/n - T-J, T 2 or T x 

according as the segment contains (1) both foci, (2) the attracting 
but not the empty focus, (3) no focus. And now, by combining 
the six cases considered separately, we arrive at the final result for 
motion within a single circuit : 

(1) t = T 1 when the segment of which the curved boundary is 

described contains neither focus. 

(2) t = 27 rjn - T x when the segment contains both foci. 

(3) t = T 2 when the segment contains the attracting but not 

the other focus. 

(4) t = 27 r/n - T 2 when the segment contains the second but 

not the attracting focus. 

4. If we are given merely the position of the attracting focus 
and of the two points 1, 2 on the orbit and the value of a, it is 
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easy to see that all the four cases enumerated are possible solutions, 
and that it is not possible to discriminate between them without 
further knowledge of the circumstances. For the second focus lies 
on the circle whose centre is 1 and radius 2 a — r Y \ it also lies on 
the circle whose centre is 2 and radius 2 a-r 2 . The necessary 
and sufficient condition that these circles shall intersect is clearly 

2a-r 1 + 2a-r 2 >c 
* i.e. 4 a>r l + r 2 + c 

c being the length of the chord joining 1, 2. 

If the condition (which, we may notice, entails the reality of 
e and 8 ) be fulfilled, the ellipse will be real, and there are two 
positions which the second focus may occupy, one on either side of 
the chord joining 1, 2. Hence in the one case the foci are on the 
same side, in the other case on opposite sides of the chord. In the 
former we have case (1) or (2) according to the arc chosen in 
passing between 1 and 2; in the latter we have case (3) or (4) 
according as the motion takes place via “perihelion,” or via 
“ aphelion.” The four cases are thus all possible, and the choice 
among them must be guided by additional knowledge. The four 
solutions are given by Gauss,* but he did not discuss the 
geometrical criteria which serve to distinguish them. Curiously 
enough, Callandreau himself omitted case (3), while the solutions 
given in the text-books are (1) and (3). When 4a = r 1 + r 2 + c the 
circles referred to above touch one another at a point between 
1 and 2. There is then only one'position of the second focus and 
the chord 12 passes through it. The four solutions become equal 
in pairs (c — i r), as they also do when the attracting focus lies on 
the chord +r 2 = c; 8 = o). 

5. Some difficulty arises in the calculation of t when e and 8 
are small or nearly equal. The former case occurs when a is large, 
i.e. when the ellipse approaches parabolic eccentricity ; the latter 
occurs when the chord e is small. Both these circumstances may 
concur. For elliptic motion I have previously f given a modi¬ 
fication of the formula which avoids any uncertainty in the 
computation, and obviates resort to any special tables. 

6. The formulae for parabolic motion can be deduced from those 
already obtained for the ellipse. As the second focus approaches 
an infinite distance it cannot lie within the finite segment whose 
boundary is described in the motion. Hence we are limited to 
the cases (1) and (3) of § 3 ; € has its principal value and 8 is 
positive or negative according as the angular motion about the 
centre of attraction is less or greater than 180°. Hence in the limit 

ae 2 = r Y + r 2 + c , a 8 2 = r 1 + r 2 - c 

6 fj)t = a£(c 8 — S 3 ) 

= (»i + r 2 + c)i + (?i + r 2 - c)i 
* Theoria Motus, §107. f M.N., lxvi. p. 496. 
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the upper sign being taken when the segment does not contain the 
focus, the lower when it does contain the focus. This is Euler’s 
equation, the use of which is much facilitated by Encke’s trans¬ 
formation and auxiliary tables. This transformation is itself not 
free from ambiguities, and the question has been recently discussed 
by S. Tscherny.* 

7. We consider next the case of motion along the concave 
branch of an hyperbola. Any point on one branch of an hyper¬ 
bola referred to its centre as origin can be represented by 

x = a cosh u , y = a(e 2 - i )* sinh u 

The focus for which the branch is concave is (ae , o). Hence 

r 2 = a 2 {(e - cosh u) 2 + (e 2 - i)(cosh 2 ?« - i)} 
r = a{e cosh u — i) 

and 

c 2 = a 2 (cosh u 2 - cosh uJ 2 + a 2 (e 2 - i)(sinh u 2 - sinh w x ) 2 
= 4a 2 sinh 2 %(u 2 - Mj){ - 1 + e 2 cosh 2 \(u 2 + u^)} 

= 4a 2 sinh 2 a sinh? f 3 

where 

a = i( u 2 " u i), cosh = e cosh \(u 2 4- u Y ) 

We may suppose u 2 >u 1 and therefore a positive; also of the 
two possible values of /•?, the positive can be chosen. Again 

r 2 + r i — 2a(cQsh a cosh /3 - 1) 

If then % € = y3 + a, 8 = 0-a 

( r i + r 2 + c)/4« = w x 2 = sinh 2 Je, (r x + r 2 - c)/4a e m 2 2 = sinh 2 

Now 

= J r 2 d0 = - ydx) 

= 2 ~ i) 4 J{(cosh ^ - e) cosh u - sinh 2 u)du 

— a 2 (e 2 — i )H( w i ~ w 2 ) - e(sinh - sinh w 2 )} 

= & 2 (e 2 - i) 4 ( - 2 a + 2 sinh a cosh / 3 ) 

— ~ i)*{ - (e - 8 ) + (sinh c - sinh 8)} 

where 

A 2 = /xa(e 2 - 1) 

yj being the attracting force at unit distance. Hence 
y}a~% t = sinh € - sinh 8 - (c - 8) 
where e and 8 are given by 

sinh ^e= ±m l , sinh 18 = ±m 2 

or 

exp. (Je) = ±m 1 + (m 2 +1) 4 , exp. (£8) = ± m 2 + (m 2 + i) 4 
* A.N., No. 4280. 
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and it is necessary to decide between these pairs of positive, 

reciprocal values. Now e is positive because a and ft are both 

positive. Hence exp. e)>i and +m 2 must be chosen. Again, 
the equation of the chord gives for the intercept on the axis 
beyond the focus 

x — ae~a sech \{u 2 + w 1 ){cosh ^(u 2 - u^) - e cosh \{u 2 + 

= - 2a sech \(u 2 + u x ) sinh Jc sinh |8 

Hence, according as the chord intersects the axis beyond the focus 
or not, i,e. according as the angular motion is greater or less than 

180 0 , 8 is negative or positive, that is, - ra 2 or + m 2 must be 

chosen. Accordingly 

fxha~H= + 2 m 1 (m 1 2 + i) 4 + 2 m 2 (m 2 2 + i) 4 

- 2 Log {(m 2 + i) 4 + mj ± 2 Log{(m 2 2 + 1) + m 2 ) 

where Log denotes natural logarithm, and the upper or lower sign 
is to be taken according as the angular motion about the attracting 
focus is less or greater than 180°. 

8. With any given (positive) value of a it is always possible to 
draw two hyperbolas such that the points 1 and 2 lie on the 
branch which is coucave to the attracting centre. For the posi¬ 
tion of the second focus is determined by the intersections of two 
circles whose centres are at 1 and 2, and whose radii are 2 a + r x 
and 2 a + r 2 respectively. The intersections are always real, for 
the necessary condition 

2 a + r 1 + 20H- r 2 >c 

is always fulfilled. They lie on opposite sides of the chord 12, 
and hence one focus must be on the same side as the attracting 
focus, the other on the opposite side of the chord. These facts 
correspond with the ambiguity of sign in the solution, which can 
only be removed, therefore, by independent knowledge of the 
geometrical circumstances of the motion. 

9. To obtain a form more suitable for calculation, we put 

m x 2 = sinh 2 = tan 2 , ra 2 2 = sinh 2 J8 = tan 2 fx 2 

where fi x and yu, 2 are to be taken in the first quadrant. Then we 
have 

tan /Xj sec / u. 1 + tan /x 2 sec /x 2 
- Log tan (45° + J j) ± Log tan (45 0 + J/* 2 ) 

In what follows only the upper signs will be retained. We 
can write 

/j)a~H = sinh c - sinh 8-2 sinh J(c - 8) 

+ 2{sinh - 8) - - 8)} 

= 4 sinh - 8) sinh 2 + 8) + J sinh 3 - 8).A 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ at University of Leeds on July 6, 2015 




1909MNRAS..69..181P 


Mr. H. G. Plummer , 


lxix. 3, 


186 
where 

A = 6{sinh £(c - 8) - |(c - 8)}/sinh 3 J( € - 8) 

and evidently tends to unity as (e-8) becomes smaller. Now 

c^a = - m 2 2 = sinh ^(c - 8) sinh 4- 8) 

and it can be found that, as a sufficient approximation when (e - 8) 
is small, 

log A = 3*6 log sech J(e - 8) 

Hence 

fj)t = ate tanb + 8) + ¥ 1 T a“$c 3 sinh -3 + 8) {sech - 8)} 3 ' 6 


and we easily find that 

tanh -j(e + 8) = sm J(/ij + /a 2 )/cos - /a 2 ) 

sinh ^(c + 8) = 2 sin + /a 2 ) cos ^(/a 2 — /a 2 ) sec /a x sec /a 2 

sech £(e - 8) = (cos cos /a 2 )* sec + /a 2 ) 

The equation has thus been reduced to a form which can 
present no difficulty to the computer when the angles /a x and /a 2 
become small or nearly equal, as they do when the hyperbola is 
nearly parabolic or the chord is small. 

io. .It has been found useful to employ the approximation 

A=6{sinh J(c - 8) - £(e - 8)}/ sinh 3 i(e - 8) 

= {sech - 8)} 3 ' 6 


and it will be well to examine the nature of this approximation. 
We have 

Log (sinh x - x) = Log (-i sin ix - x) 


= Log 


X 3 X 5 


3 ! + 5 ! + 


= 3 Log x - Log 6 + S] 


where i 2 = - i and * 


S 1 = —- 






But 


20 16800 756000 


x 2 xt 


+ . . . 


Log sinh x = Log x + j|-- 


a* 


6 180 2835 

Hence writing x for J(e - 8) we find 


Log A = - -^-£ 2 + — ^ ^ 


20 5600 84000 

Gf. A. A., No. 3883, p. 300 
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Now 


so that 


3*6 


T , # 2 x 4 x 6 

Log cosh x = -1-- ... 

2 12 45 

Log cosh \x = - -~di 4 -f -i-a; 6 • . . 

& J 20 160 800 


which shows that 


Log A — v6 Log sech hx — - — — a 4 4 -—— x° . . . 

1400 5250 

The residual terms only begin to be significant even in 7-figure 
logarithms when x is about o*i, i.e. when c — 8 is 0*2 ; and when 
the difference between € and 8 is so great there is no difficulty in 
the direct calculation of t. And it must be noted that the term 
of which A is a factor is itself small in comparison with the other 
term, and therefore A is not required with very great accuracy. 

11. The method of § 7 can now be applied to the case of 
motion on the convex branch. The focus for which the same 
branch of the hyperbola is convex is at ( - ae , o). In this case 

r = a(e cosh u+i) 


and, as before, 


c = 2a sinh a sinh /3 


where 


a = %(u 2 - u x ) , cosh f3 = e cosh ^(u 2 + w x ) 
and a and f 3 are both positive. If then 
€ = /3 + a,8 = /3 — a 


since 


r 2 + 1 \ = 2a(cosh a cosh f3 + 1) 
it follows that 


m^={r^ + r 2 + c)/ 4a = cosh 2 , m 2 1 =(i\ + r 2 - c)/4« = cosh 2 ^8 
Now 

ht — a 2 (e 2 - 1 )lj{(cosh u 4- e) cosh u - sinh 2 u}du 

= a 2 (e 2 - i)H( w 2 “ u i) + e ( 8 ^ u 2 “ u i)} 

= a 2 (e 2 - i)*{e - 8 + sinh e - sinh 8} 

where 

h 2 = fxa(e 2 - 1) 

jj, being the repulsive force at unit distance. Hence 
fJa-H = e - 8 + sinh c - sinh 8 
where c and 8 are given by 

cosh cosh £8 = + m 2 
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These are necessarily positive quantities, but two solutions are 
again possible in each case for 

exp. (Je) = m 1 ± (m 2 2 - i) 4 , exp. (£8) = m 2 ± ( m 2 2 - i)* 

these being pairs of positive, reciprocal values. As before, € is 
positive, and therefore exp. (^e)>i, and the upper sign must be 
taken in the first ambiguity. Also the intercept on the axis beyond 
the focus for which the branch is concave, in this case the empty 
focus, is as before 

x - ae = -2a sech \ (u 2 + % 2 ) sinh sinh J8 

which shows that 8 is positive, i.e . exp. (JS)>i, if x-ae is 
negative, that is, if the angular motion about the empty focus is 
less than 180 0 . Hence 

fj}a~k — 2 Log {m x 2 - 1 y + m Y } + 2 Log {m 2 2 - i) 4 4 -m 2 } 

4 - 2m 1 (m 1 2 - i) 4 + 2 m 2 (m 2 2 - i) 4 

where the upper or lower sign is to be taken according as the 
angular motion about the empty focus is less or greater than 180°. 

12. The second focus of the hyperbola in this case is one of the 
intersections of the circles whose centres are at 1 and 2, and whose 
radii are r x — 2a and r 2 - 2a respectively. In order that these 
circles may intersect it is necessary that 

7 \ - 2a -t-r 2 - 2 a>c 
i.e. 4 a<r 1 +r 2 -c 

If this condition, which entails the reality of e and 8, be fulfilled, 
two hyperbolas can be drawn so that the two points lie on the 
branch convex to the focus of repulsion. As in other cases, the 
two positions which the second focus may occupy are on opposite 
sides of the chord, and therefore in the one case it is on the same 
side as the first focus, in the other case on tfie opposite side. The 
ambiguity in the solution is inherent in the data. 

13. To reduce the last equation to a form more suitable for 
calculation we put 

m 1 2 - 1 = tan 2 />q , m 2 2 - 1 = tan 2 fx 2 

where and /x 2 are to be taken in the first quadrant. Then we 
have 

= tan {jl x sec + tan p 2 sec p, 2 
+ Log tan (45° + J/Xj) + Log tan (45° + §,u 2 ) 

The lower signs give a sum of four positive terms which 
cannot present much difficulty to calculate, and the case is not 
of practical importance. Retaining only the upper signs we find 
since 

c/2a = m 2 2 - m 2 2 = sin (/q - fx 2 ) sin (/^ + /x 2 ) sec 2 sec 2 /x 2 
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tan fa sec fx ± - tan /a 2 sec /x 2 = cl 4a. (i + sin^ sin /*, 2 ) 

~ sin 4- fc 2 ) cos J(/x, x - fx 2 ) 

and 

Log tan (45° + J/Aj) - Log tan (45 0 4- £/x, 2 ) 

= Log{cos J(/Aj + p 2 ) + sin |(/a x -^ 2 )} - Log{cos 4- /x 2 ) 

-sin 

= Log (1 4- #) - Log (1 - a:) 
where 


a; = sin - /a 2 )/2 cos + /* 2 ) cos £0^ - /x 2 ) 

= c/4a. cos 2 /*! cos 2 /x 2 /sin (p x 4- /x, 2 ) cos + y 2 ) cos £0^ - y 2 ) 

Thus a; can be calculated accurately however small and fx 2 
may be. Both parts of the expression for t have therefore been 
reduced to a form free from difficulty. 

14. If the value of x is very small, as it may be for a short 
arc, an approximation to the value of Log (1 4- x) — Log (1 — x) can 
be found thus. We have 


Log (1 + x) - Log (1 — x) = 2\ x + 4- -x 5 + - x 7 4- . 

\ 3 5 7 

But 

X 2 X 7 

tan x = x + -x s + — a: 5 4- — — x 7 4- 
3 3i5 

Hence 


2 8 

Log (1 +x) - Log (1 —x) - 2 tan x~ — x 5 4- - x 7 4 - 

I 5 45 


The approximation 2 tan (x. 206265") holds therefore with fair 
accuracy if ^<o*i; and from o*i onwards the function can be 
calculated directly under the form 


Log (1 +x) - Log (1 - x) = Log tan (45 0 + tan 1 x) 

Note that throughout Log has been used to denote natural 
logarithm, and log to denote common logarithm, so that 

1^ = 2-3025851 log = [o'3622i57] log 


University Observatory , Oxford: 
1909 January 7. 


Note added 1909 January 19. 

The approximation given in § 14 is of some importance, and 
may be useful for other purposes, e.g. in calculating the Guder- 
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mannian function for small values of the argument. It may be 
well, therefore, to examine it a little more closely. We have 

Log {Log (1 + x) - Log (1 - x)} =Log 2\x + — + — j 


= Log 2x + R 


where 


( r»2 

I+L+L+ .... 

3 5 

= lx^+— x 4 + 1 x s + . . . . 

3 9 ° 2 835 5 6 7 °° 


On the other hand, 


Log tan x = Log x + -x 2 + 1 -x 4 + x 6 + ■ 12 — x 8 + ... . 


9 ° 


2835 18900 


Hence 


Log{Log (1 +x) - Log (1 - x)} — Log (2 tan x) 


= — X 4 + — x 6 + -3 x s + .... 

15 15 5670 

= 1 5 ( 1 -# 2 ), very nearly. 

Changing to common logarithms we have then 

log{Log (1 +x) - Log (1 - x)} = log (2 tan x) 

+ [8*4616930 ]x*/(i-x 2 ) 

and the following little table gives the correction additive to log 
(2 tan x) in the seventh decimal place:— 


X. 

0*04 

Corr. 

OOI 

X. 

0*12 

Corr. 

061 

X. 

0*16 

Corr. 

195 

*05 

002 

•13 

084 

•165 

221 

*o6 

004 

•135 

098 

•17 

249 

•07 

007 

•14 

113 

•i 75 

280 

•08 

012 

•145 

131 

•18 

314 

•09 

019 

•15 

150 

•185 

351 

•10 

029 

•155 

171 

•19 

391 

•11 

043 

•16 

195 

•195 

435 

•12 

061 



•20 

483 


Thus, for example, when x — 0’2, which is the circular measure of 
11° 27' 32 / ' , 96, we have log (Log i* 2 /o* 8 ) = 9*6079536 and log 
(2 tan x) = 9*6079053, the difference coinciding exactly with the 
last entry in the table. 
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A Method for the Reduction of Comet Photographs» 

By H. C. Plummer, M.A. 

1. The magnificent results which have lately been obtained in 
comet photography suggest some reflections on the method by 
which the plates may be reduced. The method of Bessel,* 
which was designed for the treatment of Halley’s comet in 1835, 
and has been followed more or less faithfully by all subsequent 
writers, deals with observed position-angles and distances with 
respect to the nucleus. If the observations are not made in this 
form, they must first be reduced to it. Of course position-angles 
and distances can be very easily deduced from the photographic 
plate, f But it has frequently been found | that there is an 
advantage in the use of rectangular coordinates over spherical 
geometry in dealing with photographs. There is a corresponding 
change in the method of calculation from the general use of 
logarithms to the calculating machine which appeals strongly to 
some computers. It may, therefore, be worth while to investigate 
the formulae which seem to be the natural expression of this point 
of view. Whether the emanations from a comet remain strictly in 
the plane of the orbit is open to doubt, but provisionally it is 
necessary to make this assumption, and it forms the basis of what 
follows. 

2. It will be convenient to give precepts for certain preliminary 
operations. Let us form the arrays 

- sin A cos A o 

- sin D cos A — sin D sin A cos D (1) 

cos D cos A cos D sin A sin D 


where A, D are the R.A. and declination of the plate-centre, and 


X 

a sin (a + u) 
a cos(a + u) 


Y 

b sin Q 3 + u) 
b cos (/3 + u ) 


Z 

c sin (y -I- u) 
c cos (y + u) 


<*> 


where X, Y, Z are the rectangular coordinates of the Sun referred 
to the equatorial system of coordinates, a, b, c, a, /?, y the Gaussian 
constants for the orbit of the comet, and u an angle measured from 
the ascending node in the direction of the comet’s motion: the 
actual value to be assigned to u can be left for definition later on. 
Next we form the array 



n 

m 2 

n \ 

m B 

n. 


( 3 ) 


* A.N., No. 300. t Of. M*N., lxix. p. 100. 

t Cf. M.N., lx. pp. 176 and 201 ; lxii. p. 29. 
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